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Abstract 

From the perspective of the network theory, the present work illustrates how the 
parametric intrinsic geometric description exhibits an exact set of pair correction functions 
and global correlation volume with and without the inclusion of the imaginary power 
flow. The Gaussian fluctuations about the equilibrium basis accomplish a well-defined, 
non-degenerate, curved regular intrinsic Riemannian surfaces for the purely real and the 
purely imaginary power flows and their linear combinations. An explicit computation 
demonstrates that the underlying real and imaginary power correlations involve ordinary 
summations of the power factors, with and without their joint effects. Novel aspect of the 
intrinsic geometry constitutes a stable design for the power systems. 
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1 Introduction 



Planning issues towards a stable power supply, viz., transmission and distribution systems 
have been in wide application of the power system since a decade [T]. In the transmission 
theory, this shows that the power flow varies as a function of the power factor of the network. 
In order to regulate the voltage of the network, the power factor (p is thus related to the 
impedance angle a > ~ — (p. Interestingly, the determination of a can be accomplished by 
tuning the network parameters, viz., resistance (r), inductance (L) and capacitance (C) [2j[3]. 
An optimal choice of a improves the efficiency of power networks, and thus an appropriate 
design of the network parameters. For a given transmission line(s), we provide a proper network 
planning for the operations with a minimum reactive power requirement. From the outset of 
the intrinsic geometry, the present paper determines the required unit of the power flow under 
the fluctuations of the power factors. 

Up to now, most of the network planning and design is based on the power flow equa- 
tions jU[S] , and thus the network characterizations are linearly afforded by a set of chronological 
data analysis, heuristics methods, parametric estimations and optimization techniques |U[7]. 
In this paper, this motivates us to define non-linear criteria for the electrical networks. Subse- 
quently, we take an account of the fact that the set of voltages at all buses attains an equilibrium 
configuration, therefore the extremization of the power flow in the network can be determined 
in terms of a. The selection of the network parameters determined by the optimization tech- 
niques could be unreliable and thus the cause of bottle-necking. This shows an urgent need for 
compensator (s). For a safe operation and optimal power flow, our analysis provides stability 
criteria for the non-linearly reliable comportment of the electrical network. In this concern, our 
method provides a high degree of compensation strategy to reduce the fluctuation effects of the 
network parameters, viz., r, L and C. 

To be specific, we focus our attention on the power networks and determine the required 
set of voltage stability criteria, viz., the selection of power factors, network planning and com- 
pensation strategies. Further, the intrinsic geometry offers an effective determination of the 
power system characterization, power factor corrections and voltage regulation. It is worth 
mentioning that the voltage level at the buses could be tuned in the equilibrium. In the reverse 
engineering, the proposed solution keeps the power system in a stable voltage range, under the 
fluctuations of the network configuration. In this way, we provide an efficient power system 
characterization, which is non-linearly stable over the fluctuations of the network phases. 

In this concern, the intrinsic geometry has been important in the configurations involving 
black holes in string theory [9TTT5] and M-theory [T6HT9] . possessing a set of rich stability struc- 
tures [2"0H2"5] . Therbay, there have been several investigations about the equilibrium perspective 
in black holes, explicating the nature of pair correlations and the associated stability of the so- 
lutions [31]. Besides several general notions analyzed in condensed matter physics [291433] . we 
consider specific electrical network(s). As mentioned above, we analyze the parametric pair 
correlation functions and their correlation relations about the equilibrium configuration. We 
find that the intrinsic geometric consideration entails an intriguing feature of the underlying 
fluctuations, which are defined in terms of the network parameters. 

Given a definite covariant intrinsic geometric description of the network configuration, we 
expose (i) the conditions for the stability, (ii) properties of the parametric correlation functions 
and (iii) scaling relations in terms of the parameters of the network. In this analysis, we enlist 
the complete set of non-trivial parametric correlation functions of the electrical networks. This 
anticipation follows from similar considerations in the black hole solutions in general relativity 
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[?, ?,E3ES], attractor black holes [3HH52] and Legendre transformed finite parameter chemical 
configurations [53j[5l], quantum field theory and associated Hot QCD backgrounds [55],[56], as 
well in the stability of Quarkonia states [57]. Thus, the intrinsic geometry plays an important 
role in the study of power networks and their stable design applications. 

In the context of power flow equations [HEj , the usefulness of the present investigation may 
apparently seem to be limited, however, we illustrate that our preposition could in principle be 
generically applied to all electrical networks, in order to achieve the best understanding of the 
phenomenon and importance of the controlled power flow and network analysis. The intrinsic 
geometric consideration is capable to provide strategic planning criteria for the effective use of 
power system and voltage stability. We show that this notion follows from the standard laws of 
the electrical circuits [5]. Subsequently, for an additional component, the criteria of the voltage 
stability can be used for an optimal selection of the network parameters. 

2 Phases of the Electrical Network 

In this section, we set up the formulation of the problem and outline the notion of the intrinsic 
Riemannian geometry. The exploration of the power flow equations gives the relation of the 
power flow with network phases, and thus the consideration for the analysis. Most efficiently, 
our intrinsic geometric model is designed to provide the critical values of the phases. It is worth 
mentioning that the present method is important towards the determination of the parameters, 
and thus the unstable mode of the network. 

2.1 Hypothesis of the Power Flow 

For the optimization of the power flow, we share the hypothesis of the prior solutions and use 
the load flow equations jH[S] in order to solve the mentioned issues. The associated power 
conservation equations with the real (resistive) and imaginary (reactive) branch parameters are 
given by 

Pi = ^2\Vi\\Vj\\Yij\(GijCos(a i: j + Sj - Si) + Bijsin(a,ij + Sj - Si)) 

(1) 

Qi = ^2\Vi\\Vj\\Yij\(Gijsin(a,ij + Sj - Si) - BijCos(aij + Sj - Si)) (2) 

In the above equations, the phases are defined as 

X r . . — Xn.. 

tan{aij) = — ^ ^ (3) 

The inverse set of the impedances Zij and voltage angles 5j are 



(ry+PfLy -jX Cij ) 
Sj = A (4) 

m 

For the purpose of subsequent analysis, let us consider arbitrary i th -hv& such that the 
steady state condition is realized as \Vi\ = 1. This makes the underlying configuration reach an 
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equilibrium. Thus, the respective cases of the present interest reduce to the standard network 
considerations. For the purpose of the future analysis, we consider that a loss-less line is defined 
by ciij = 90°, where +90° represents the ideal case. In this case, it turns out that the network 
is purely inductive with r = 0. Notice that a realistic network would never reach the limit of 
the zero resistance. In the generic situations, the value of the varies from +90° to —90°. 
Further, the phase = —90° is also not feasible in the real situation, as the network possesses 
a finite capacitance. Thus, the phases for the inductor and capacitor circuits can be defined as 

— 1 / ^ L a n 

a (l)ij = * an ( )l a ij = — 90 

Xc 

a {2)lj = tan' 1 (— 2-), Oij = 90° (5) 



For a general consideration of the network fluctuations, we have non-zero values for the 
network parameters, viz., r, L and C, and thus the general phase angle a^)ij is defined as 

o (3)ij - = tan-\ ^) (6) 



We take an account of the fact that the efficiency of the power flow on a transmission line 
is analyzed by the phases of the impedance pertaining to the transmission lines. As per the 
consideration of the next subsection, our method offers a non-linear characterization for the 
generic component of a realistic network, which we suppose neither a purely inductive nor 
capacitive component. 



2.2 Hypothesis of the Intrinsic Geometry 

In this subsection, we recall the motivation for intrinsic geometric analysis and set up the 
notations for the subsequent computations. Following the notations of the previous subsection, 
a given network can reach a local equilibrium, if we can fix one of the phases of the power 
network. The logic simply follows from the fact that the sum of the three angles of the trigone 
is a constant. To be specific, let us illustrate the consideration of the intrinsic geometry for 
the case of two parameter configurations. To be concrete, let the parameters be {01,02} and 
let S (01,02) be a smooth function of the network (real, imaginary ) phases as defined in the 
Eqns.([U E]) or any of their real combinations. For a given S (0,1,0,2), the components of the 
correlation functions are described as the Hessian matrix Hess(S(ai, 02)) of the generalized 
power function under the flow of the parameters. Following this consideration, the components 
of the intrinsic metric tensor are given by 

d 2 S d 2 S d 2 S 

9aiai ~ doT 9aia2 ~ doTdo? 9a2a2 ~ daj { ' 

The components of the intrinsic metric tensor are associated to the respective pair correlation 
functions of the concerned power flow. It is worth mentioning that the co-ordinates of the 
underlying power factor lie on the surface of the parameters, which in the statistical sense, 
gives the origin of the fluctuations in the network. This is because the components of the 
metric tensor comprise the Gaussian fluctuations of the network power, which is a function 
of the parameters of the power configuration. For a given network, the local stability of the 
underlying system requires both the principle components to be positive. In this concern, the 
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diagonal components of the metric tensor, {g ai(H | i S 1,2} signify the heat capacities of the 
system, and thus they are required to remain positive definite quantities 

9aiai > 0, i= 1,2 (8) 

From the perspective of intrinsic geometry, the stability properties of the network flows can 
thus be divulged from the positivity of the determinant of the metric tensor. For the Gaussian 
fluctuations of the two charge equilibrium power configurations, the existence of a positive 
definite volume form on the power surface imposes such a stability condition. Specifically, a 
power supplying configuration is said to be stable if the determinant of the tensor 

||<?|| — S aiai S a2a2 — S aia2 (9) 

remains positive. For the two parameters networks, the geometric quantities corresponding to 
the chosen power elucidate the typical features of the Gaussian fluctuations about an ensem- 
ble of equilibrium states. Subsequently, the intrinsic scalar curvature, as a global invariant, 
accompanies the information of the correlation volume of the underlying power fluctuations. 
Explicitly, the scalar curvature R takes the following form: 

R = ^ {^a 1 a 1 S a2 a 2 &a 1 a 2 ) (Sa 2 a 2 S ai a iai S aia2a2 



, q q q , q q q 

~ k -'aia2 kJ aiaia 2 kJ a\a2a 2 < ^ajaj ^a-^ai 0,2 ^ 0,20.20.2 

'aiai 010,20,2 0202'-' a\a\a 2 



Saia 2 ^aia 1 a 1 S a2a2a2 S aiai S a a a S a2a2 S a a a ) (10) 



Notice that the zero scalar curvature indicates that the power of the network fluctuates inde- 
pendently of the phases, while a divergent scalar curvature signifies a sort of phase transition, 
indicating an ensemble of highly correlated pixels of information on the power surface. In the 
case of black hole physics, Ruppeiner has interpreted the assumption "that all the statistical 
degrees of freedom of a black hole live on the black hole event horizon" as an indication that 
the state-space scalar curvature signifies the average number of correlated Planck areas on the 
event horizon of the black hole [29]. For the case of the two parameter systems, the above 
analysis of the surface shows that the scalar curvature and curvature tensor are related by 

2 

R(ai, CI2) = Y^Ra 1 a 2 a 1 a 2 (H) 

The scalar curvature thus defined offers the nature of the long range global correlation and 
underlying phase transitions originating from the power flow. In this sense, we anticipate that 
an ensemble of signals corresponding to the network are statistically interacting, if the underly- 
ing power configuration has a non-zero scalar curvature. Incrementally, we may notice further 
that the configurations under present consideration are allowed to be effectively attractive or 
repulsive, and weakly interacting, in general. The intrinsic geometric analysis further provides 
a set of physical indications encoded in the geometrically invariant quantities, e.g., scalar cur- 
vature and other geometrically non-trivial objects. For the electrical network, the underlying 
analysis would involve an ensemble or subensemble of the equilibrium configuration forming a 
statistical basis about the Gaussian distribution. With this brief introduction, we shall now 
proceed to systematically analyze the underlying stability structures of the network fluctuations 
in the real, imaginary power flows and their joint effects on the network. 
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3 Real Power Flow 

Let us first describe the intrinsic stability of the electrical network with a given power factor. 
Following the Eqn.Q, the power defined with a set of desired corrections over the network 
power factors, chosen as the network variables ai, 02 for the present analysis, is given by 

P( Ql ,a 2 ) := — — - V \ (12) 
Ko(l + (tan(ai) — tan(a2)y) 

The components of the correlation functions are described as the Hessian matrix Hess(P(ai, 02)) 
of the concerned power under the tuning response function. Following the Eqn. (ll2p . the com- 
ponents of the metric tensor are 



Qa\a\ 



2c\V 2 nfj 
Ro 

^ 2 c\V 2 nf 2 
Ro rf 2 
2cf V 2 n£ 



( 13 ) 

' 122 



In this framework, we find that the geometric nature of the parametric pair correlations offers 
the notion of fluctuating networks. Thus, the fluctuating parameters may be easily divulged in 
terms of the intrinsic parameters of the underlying network configurations. For a given network, 
it is evident that the principle components of the metric tensor signify self pair correlations, 
which are positive definite functions over a range of the parameters. In order to simplify the 
subsequent notations, let us define the cosine and sine variables as 

cos(ai) := Ci,i = 1,2 

sin(ai) := Sj,i = 1,2 (14) 



We have arrived at the conclusion that the numerator of the local pair correlation functions 
are expressed as the following trigonometric polynomials: 



nf 2 := -6cfc 2 + 6sic|s 2 cf - c| + 6cfc| 

— 2s\S2c\ + 3cfc2 — c\c2 
n ?2 := 6siC2S2Ci — 3cf + 7c\<% — 3c| 

"22 : = ~ c l c 2 - c l - 2S1CJS2 + 6c2«lS2cf + &c\c2 

-Gcxci^ + Scxc 2 . (15) 



We further notice a similar conclusion for the denominator of the local pair correlation func- 
tions. In general, we find, for the real power flow pair correlations, that the denominator of the 
local pair correlation functions are 



r « — - c f - 15cf 4 + A2c\c\ - 27cfc| + 15c? csj - 27c?c| 

+15cf c| + 6sic 2 S2cf + 20sic|s2cf - 15cf c| 

+6sic|s 2 ci - c\ + 13c? - 20sic^s 2 c? 

-20sic|s2C? + 14c?C2Sis 2 
rf 2 := -c? - 15cf c| + 42cf c\ - 21c\c% + 15c? csj - 27c? c| 

+15cf c^ + 6sic 2 s 2 c? + 20sic|s 2 c? - 15cf c| 

+6sic|s 2 ci - c\ + 13c? cj] - 20sic^s 2 c? 

— 20sic|s2C? + 14c?c|siS2 
r 22 : = ~ c i ~ 15c l °2 + 42c i c 2 - 27cfc§ + 15c?c^ - 27c? c^ 

+15cf c| + 6sic 2 s 2 c? + 20sic|s 2 c? - 15cfc| 

+6sic|s 2 ci - cf + 13c? c§ - 20sic|s 2 c? 

-20sicjjs 2 Ci + 14c?cjjsis 2 (16) 



5 



It is worth mentioning that the real network is well-behaved for the generic values of the 
parameters. Over the domain of the {01,02}, we observe that the Gaussian fluctuations form 
a set of stable correlations, if the determinant of the metric tensor 

8V 4 cfc? n R 

9= J 1 17 
R r g 

remains a positive function on the power factor surface (M 2 (i?),g). In terms of the trigono- 
metric polynomial, we obtain that the numerator of the determinant of the metric tensor can 
be expressed as 

R. - 10„. „_„4„4 _ A „3„3 _|_4 C 3 C2 _ C 5 C2 _ ClC 5 



a 



+4cicf — c\s\S2 — cfsiS2C2 — cf sis 2 c? 
~4 s i s 2 — 7c? cf — 6c\c^s\S2 — 7e?cf 

+12c?cf (18) 

As per the expectation, the denominator of the determinant of the metric tensor turns out to 
be given by the following trigonometric polynomial: 

rf : = -21044 - 210c?cf + 910c?cf - c?° - cf + 121c 2 °c?° 

+122efc?sis 2 + 252cf c|sis 2 - 584si4s 2 c? + 332sic 2 ! s 2 c? 
+808sicfs 2 c? + 120sic|s 2 ci - 120cfc?sis 2 + 10sis 2 c 2 c? 
+45c\°c\ - 2b0c\°c\ + 490cf°c? - 405^00? + 45c|c?° 
-120siC2S 2 cf + 10sic 2 s 2 ci - 584sic 2 s 2 Ci - 344sic§s 2 Ci 
-344si4s 2 cf + 332cf c?sis 2 + 120si4s 2 c? - 250cfc?° 
+460cf c? - 45c? cf + 490c|cJ° - 1190c|cf + 1180c|cf 

-405cf cjo - 1190c|cf - 45cfc? + 460cf c? (19) 

Herewith, the behavior of the determinant of the metric tensor shows that such a real power flow 
becomes unstable for the specific values of the parameters. For generic Rq and V, the nature 
of the determinant of the metric tensor is depicted in the Eqn. (fTT|) . It is worth mentioning 
further that the electric networks become unstable in the limit of vanishing {ai, 02}. 

In order to explain the nature of transformation of the {a\, a 2 } forming the intrinsic surface, 
let us explore the functional behavior of the associated scalar curvature. Our computation shows 
that the scalar curvature reduces to the following form: 

1 

R = 4R c? C2 V 2 rj (20) 

The numerator and denominator of the scalar curvature take the following trigonometric 
polynomial expression: 



r4 : = 



-42c? cf - 42c? 4 + 444c? cf + 3c?° + 3c^° + 702c^°cf 
+702c 2 '4s 1 S2 + 84c?cf sis 2 - 280sic r ! s 2 c? + 356sicf s 2 c? 
-fl080sicj,s 2 c? - 24sicfs 2 c? + 24cfc?sis 2 - 18sis 2 c 2 e? 
-384°cf - 704°cf + 752c 2 °c? - 1349cf°e? - 38cfc?0 
■f24sic|s 2 cf - 18sic|s 2 ei - 280sicfs2e? - lOOOsicfsijc? 
-lOOOsicf^e? + 356cf c?sis 2 - 24sic^s 2 c? - 70ef c}° 
+7244 + 39c? c| + 752cfe?0 - 1082cfc? + 2288cfc? 
-1349cfc}° - 1082cfc? + 39cfe? + 72cfc? 

w curvature can be expressed as the trierm 



(21) 

The denominator of the scalar curvature can be expressed as the trigonometric polynomial 

rjj := -c? + 6sic 2 s 2 c? + 20sicf s 2 c? + 6sicfs 2 ci + 8sicfs 2 c? 



— CJ + OSlC 2 S 2 CJ +■ ZUS1C2S2CI + OSlC2« 2 Cl + »SlC2« 2 C 1 

+8sicf s 2 c? - 15c? cf + 6c? 4 + 53c? cf + 4c? c| + 53c? cf 
+4c?cf - 15c?cf - 48e?cf + c? + cf - 70c? cf s lS2 - cf 
-12sicfs 2 c? + 72sicfs 2 c? - 12siefs 2 e? - 12sicf s 2 c? 
-12sic^s 2 c? - llcf e? + 2c? cf - 48cf c? + 72cf c? 
-48cf c? + 2cf c? - llcf c? 



(22) 
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Figure 1: The determinant of the metric tensor plotted as the function of the power factors 
ai, a 2 , describing the real power fluctuations in electrical networks. 




Figure 2: The determinant of the metric tensor plotted as the function of the equal power 
factor a := ai — a 2 , describing the real power fluctuations in electrical networks. 
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Figure 3: The curvature scalar plotted as a function of the power factors 01,02, describing the 
real power fluctuations in electrical networks. 
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Figure 4: The curvature scalar plotted as a function of the equal power factor a 
describing the real fluctuations in electrical networks. 
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We find that a typical real power network is globally correlated over all the generic Gaussian 
fluctuations of the parameters {011,02}, unless = 0. For the above real power networks, we 
observe that the scalar curvature diverges in the limit = 0, showing a signature of the global 
instability on the {01,02} surface. Thus, the intrinsic geometric analysis shows that a real 
power network is interacting and locally stable over the surface of fluctuation, if the network 
parameters {ai, 02} are properly chosen. 

For the choices V = 1 and Rq = 1, the Fig. ([I]) shows the determinant of the metric tensor. 
These plots explicate the nature of the stability in the real power flow in the power networks. 
The corresponding plot for the scalar curvature is depicted in the Fig.(|3]). This plot shows 
the global nature of real power flow in the electrical network under the effects of Gaussian 
fluctuations of the parameters. 

For the equal phases, viz., a\ = a and 02 = a, the surface plots of the determinant of the 
metric tensor and scalar curvature are respectively shown in the Figs.© and (j4j). We observe 
that the stability of the real power networks exists in certain bands. This follows from the fact 
that the instability is present only for a set of specific equal values of the parameters. For a 
limiting equal phases network, the limiting scalar curvature interestingly simplifies to the shown 
shape. In particular, the corresponding peaks in the Figs.(j3]) and fll]) of the scalar curvature 
indicate the graphical nature of the global instability in the real power networks. Physically, 
the peaks in the curvature show the presence of non-trivial interactions in the network. 



4 Imaginary Power Flow 

In the present section, we analyze the nature of an ensemble of fluctuating electrical networks 
generated by a pair a\, 02. To focus on the most general case, we chose the variable a±, a% as a 
function of L, C and r of the chosen network. Following Eqn.((2]), when the imaginary power 

Q {aua2):= Yl (tan(a,)-tan{a 2 )) 

Ro (1 + (tan(ai) — tanian)) ) 

is allowed to fluctuate as a function of the 01,02, we may again exploit the definition of the 
Hessian function Hess(Q(ai, 02)) of the imaginary power. Herewith, the components of the 
metric tensor are given by 



9a±a 2 



2c 2 V 2 




Ro 


' 11 


2c 2 c 2 V 2 


n 'l2 


Ro 


r 1 
' 12 


2c? V 2 


n 22 


Ro 


r 1 



(24) 

In this case, the numerators of the local pair correlation functions are expressed as the following 
trigonometric polynomials: 

rafx := — 5c|cf S2 — sicf + 8c|sicf — 4c2«2cf + 3c2S2cf 

— 'Ac 2 s\C\ + 8cf S2C 2 — 7s\c\c 2 + S2C 2 + c 2 sic\ 
n 12 := c 2 s\ + 3c 1 siC2 — C{S2 — Tc^^si + 7c 2 S2Cf — 3e 2 S2Ci 
n 22 := — Sc 2 c 2 S2 — sicf + 5c 2 sic± — C2S2cf + 3C2S2C 2 

-3t|sici + 7cf s 2 c| - 8sicfc| + s 2 cl + 4,c 2 s 1 c 1 (25) 
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While, the denominators of the local pair correlation functions take the following trigonometric 
expressions: 



rfj := 1344 + 14cf c|sis 2 + 6sic|s 2 ci + \h44 

+6sic 2 s 2 c'i - 20c?c 2 sis 2 - 20cfc 2 sis 2 - 15c? e 2 
+20c?c 2 1 sis 2 - 21c\ 4 + 15c?c§ - 27c? c^ + 42cf cf 

— 1544 — c 2 — c? 

r{ 2 := 13c? cl 1 + 14cf c|sis 2 + 6sic|s 2 ci + 15c? 4 

+6sic 2 s 2 c? - 20cf c 2 sis 2 - 20c?c 2 sis 2 - 15c? e 2 
+20c?clsis 2 - 27cf 4 + 15c? c^ - 27c? c| + 42cf c^ 

— 1544 — c 2 — c? 

r| 2 := 13c? c^i + 14cf c|si« 2 + 6sic|s 2 ci + 15c?c^ 

+6sic 2 s 2 cf - 20c?c 2 sis 2 - 20c?c|sis 2 - 15c? ef 
+20c?c 2 1 sis 2 - 21c\4 + 15c? 4 ~ 2744 + 42cf cf 

-1544 -4-4 (26) 

It follows that the pure pair correlations {^11,^22} between the parameters {01,02} remain 
positive, which is the same as for the flow of the real power. A straightforward computation 
further demonstrates the over- all nature of the parametric fluctuations. In fact, we find that 
the determinant of the metric tensor reduces to the following expression: 



44 n l 



R 2 



(27) 



In this case, the numerator of the determinant of the metric tensor is given by the following 
trigonometric polynomial: 



n\ : = 1044 + 37c?e 2 'sis 2 - 1,44 ~ 3c i c 2 ~ 39c? 4 + 28c i 4 

+ 10e?ef - 39cfc| + 3sic|s 2 ci - 20c?c|sis 2 - 20c?c§sis 2 

+3sic 2 s 2 cf + 2c?c|siS2 + 38cfcf - 4 ~ 4 (28) 

Explicitly, we find that the denominator of the determinant of the metric tensor can be 
presented as 

r'g : = 252cfc|sis 2 + 910c? 4 ~ 210cfc| - 210c? ef 

-4° ~ 4° + 122c?c^sis 2 + 121cJ°c^° + 10sic|s 2 ci 
-120s!C^s 2 c? - 120c?c|sis 2 + 332c? c^s^ + 332c?c|sis 2 
-344c?c^sis 2 - 344c?c^sis 2 + 490c? 4° ~ 405c? 4° ~ 250cJ°c 2 
+490ci°c^ + 1180c|c? - 1190c|c? - 1190c|c? - 250cf 4° ~ 405ci°cl 
-45c? 4 ~ 4544 + 460c? 4 + 460c|cf + 45c}°c 2 + 45c 2 °c? 
+120c?c|sis 2 - 584c?c|sis 2 + 120c?c^sis 2 - 584e?e 2 'sis 2 

+10sis 2 c?c 2 + 808c?c^sis 2 (29) 



It is not difficult to compute the exact expression for the scalar curvature describing the 
global parametric intrinsic correlations. In particular, we find that the scalar curvature reduces 
to the following form: 

R --\~4^~ 2 V R (30) 

In the above equation, the numerator of the scalar curvature takes the trigonometric polynomial 
expression 

n^ 1 := -297sic?cl + 2235sic?c^° + 4770sic^c| - 15596sic?c^° 

-19538sic?cl + 45624sic?c^° - 7857sieJ 1 c 2 1 - 3954s lC ?e^ 2 
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Figure 5: The determinant of the metric tensor plotted as the function of the power factors 
ai,a 2 , describing the imaginary power fluctuations in electrical networks. 
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Figure 6: The curvature scalar plotted as a function of the power factors a\, a2„ describing the 
imaginary power fluctuations in electrical networks. 



+19890sicJ4 2 - 49482sicfc^ 2 - 947s 2 cfc 2 l + 59571sic^4 2 
+121s 2 c\°4 ~ 2235s 2 c}°c 2 ' + 30498sic} 1 cf - 3653s 2 cfc| 
-lls 2 c} 4 c 2 + 15596s 2 ci°4 - 45624s 2 c}°c 2 ' + 60771s 2 c}° cl 1 
-355s 2 c} 2 c 3 > + 3954s 2 cJ 2 cl - 19890s 2 c} 2 c 2 ' + 49482s 2 c} 2 c 2 ' 
-59571s 2 c} 2 c 2 1 - 60771sicJ 1 4° + 19538s 2 c? c§ + 7857s 2 c?c 2 1 
-30498s 2 cf4 1 + 2080sic^4 4 - 8902sic}c 2 4 + 20139sic?c 2 4 
-22407sic} 1 c 2 4 + 8902s 2 c} 4 c 2 ' - 20139s 2 c} 4 c 2 l + 22407s 2 c} 4 4 1 
n^ 1 : = lls!Ci4 4 - 2080s 2 cJ 4 c| - 243sicf4 4 - 9450s 2 c 2 3 c} 4 
+61s 1 c\ 3 c 2 1 - 812sic} 3 c| + 5044sic} 3 c^ - 16769sic} 3 c| 
+30165sic} 3 4° - 27138sic|[ 3 4 2 + 9450sic} 3 c 2 4 
-5044s 2 cfc^ 3 - 61s 2 c 2 4 3 + 220s 2 c 4 c 2 ' - 121sic?c 2 ° 
+9A7s 1 c\ 1 4 + 355sicf4 2 - 45sic} 1 c 2 . + 27138c} 2 c 2 l3 .s 2 
+330s 2 c?c| + 3653sic?cf - 220sic?cf - 10s lCl cf 
+297s 2 c?c| - SZOs^lcl + 45s 2 c 2 cl 1 + 10s 2 c} 2 c 2 
-4770s 2 cf c 7 2 - 30165s 2 cJ°4 3 + 812s 2 e 4 4 3 - s lC } 3 

+16769s 2 cf4 3 + s 2 cf + 243s 2 c} 4 c 2 1 (31) 

The function appearing in the denominator of the scalar curvature can be written as the 
following polynomial: 

rjj := 42cfc| + 4 2 +cP + 5344cfc^is 2 + 10447ci 4 +4383cl 2 cf 
+100c?4 1 si«2 - 708cf4 1 sis 2 + 100c} 1 c|sis 2 + 2528c}4 1 sis 2 
-4406c?4 1 sis 2 + 2528c} 1 4sis 2 - 4406c} 1 cfsis 2 + 2812c} 1 c^ 1 sis 2 
-6cic^sis 2 - 708c\ 1 4s 1 s 2 - §c\lc 2 s 1 s 2 - 29c\c\ 2 - 29c} 2 c| 
+4383cfc 2 2 - 5813c}°c^ 2 + 2813e} 2 4 2 - 1598cf c\ 2 - 1598c} 2 c| 
-5813c} 2 c^° + 307cf c\ 2 + 307c} 2 c| - 22sic 2 's 2 cf - 22c?c|sis 2 
+368cf c 2 ) s 1 s 2 + 368c?c|sis 2 - 2132c? c\sis 2 - 2132c}c|sis 2 
+1826cf c\° - 6442cf c\° - 257c}°c| + 1826c}°c| + 2822c§c? 
-4B2c|cf - 482cf e? - 257cf4° - 6442c}°c| + 27cf c\ + 27cfc 4 

+15c}°cl + 15c 2 °c 2 - 36cic 2 sis 2 - 36c5 i c 2 sis 2 + 472c}c 2 sis 2 (32) 

Consequently, we may easily analyze the underlying conclusions for the specific considerations 
of the variable power factor of the network. As in the case of real power flow, the global nature 
of scalar curvature and associated phase transitions can be thus determined over the range of 
power factors describing the important network of interest. 

As in the previous subsection, we shall focus our attention on same electrical network and 
on the specific values V = 1 and Rq = 1. The determinant of the metric tensor shown in the 
Fig. §5$) describes the phenomenological property of the Gaussian imaginary power fluctuations. 
The scalar curvature as shown in the Fig. (JTj) depicts a couple of antisymmetrical fluctuations 
for the imaginary network power flow. 

For the equal values a\ — a and a-i = a, we notice, from the Fig. ([6]), that the system acquires 
a couple of locally chaotic fluctuations. Herewith, we find the surprising fact that the imaginary 
power flow has a different geometric nature, and the scalar curvature turns out to be zero in 
the limit of the equal phases for the imaginary power flow. It is worth mentioning that such a 
power flow is globally non-interacting. 

5 Complex Power Flow 

In order to further understand the nature of generic electrical networks, we shall now consider 
the linear combination of the real and imaginary power flows in the network. The associated 
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Figure 7: The curvature scalar plotted as a function of the power factors ai, 02, describing the 
imaginary power fluctuations in electrical networks. 



joint power flow F takes the following form: 

F( ai ,a 2 ):=Yl {l + tan{a,)-tan{a 2 )) 
Ro (1 + (tan(ai) — tan(a,2)) ) 

In order to obtain the components of the metric tensor in the power space, we employ the 
definition of the Hessian matrix as indicated previously, and see that the metric tensor has the 
following expression for the components: 



2c 2 2 V 2 ng 
Ro rg 
2c 2 qV 2 ng 2 
Ro ^ 
2c 2 V 2 n 



V> 



c 



ao.20.2 = c~ ( 34 ) 

Ro r% 2 

In the above equation, the numerators of the local pair correlation functions are expressed as 
the following trigonometric polynomials: 

™fl : = - c 2«2 +sic? + Zc\c% - 2c 2 sis 2 c? + 5clc%s 2 

+4cfc 2 s 2 - cic|si + Gc\c 2 s\s 2 — 3c 2 c 2 s 2 

+7cf4 *i - 8cf c|s 2 + 3ciclsi + 6cfc| 

-6cfc§ - 8c? cf si - efc? - c| 
n i2 := 6cjc|siS2 — 3c? c 2 — 3cf S1C2 + c?s 2 — 3cic 2 

+7e?c 2 1 + ls\c\(? x - 7s 2 c{(? 2 - c\s\ + 3c|s 2 ci 
n 22 := — c|«2 + sicf + 3cf c| + 8cf C 2 S2 + cf C2S2 

— 4cic|si + 6c?c 2 sis 2 — 3cjC2S2 + 8c?c 2 t si 

—7cfc 2 s 2 + 3cic 2 si + 6cfc 2 — c\ — cfc 2 

-5cfc|si - §c\(? x - 2cls 1 s 2 c 1 (35) 

The denominators of the local pair correlation functions reduce to the following trigonometric 
polynomials: 

rg x := -2044s lS2 - c% - cf - Tlc\c\ + 14c?c|sis 2 



13 



+13c? 4 + 6sic|s 2 ci + 15c? 4 ~ 27cj4 
-20c? c|sis 2 + 6c 2 c?siS2 + 20c?c|sis 2 + 42c? ef 
+15c?c 2 : - 15e?c| - 15c|c? 
rf 2 := -20clc?sis 2 - c| - c? - 27cf c? + 14c?clsis 2 
+13c?c| ; + 6sic!s 2 ci + 15c? 4 ~ 2744 
— 20c?c|si«2 + 6c2C?siS2 

+20c?c 2 1 siS2 + 42c? 4 + 15c?cl - 15c? 4 ~ 15cf c? 

rg := -20clc?siS2 - 4 ~ 4 ~ 27cfc? + 14c? c^U^ 

+13c?c| + 6sic|s 2 ci + 15c? 4 ~ 27cf4 

-20c?c|sis 2 + 6c 2 c?sis 2 + 20c?C2SiS2 
„4„4 , 1K „6„2 1K „4„2 1C r„4„2 



+42c{4 + 1544 - 1§44 - 1544 (36) 

The determinant of the metric tensor turns out to be a rational polynomial function in the 
{a^a^}, which in turn is given, in compact notations, as 

4V 4 c 2 c?n£ , 
9 = sf^-gr 37 

Herewith, the numerator of the determinant of the metric tensor is given by the following 
trigonometric polynomial: 

- c ■ = -2544 ~ 18c|c?sis 2 - 4 ~4 + 13c?c|siS2 



9 ■ 

— c?si + S2c| + 5S1C2S2C1 — 18c?C2SiS2 + 5c2C?siS2 
+V2c\4 - 25c? 4 + 12 4 c 2 + 5c?c 2 s 2 - 5c lC f Sl 
-10c?c|si - llc?c| - llcfc? + 10c?c|s 2 - 60sicfc? 
+60s 2 c|cf - 20c? 4 S 2 + 20e?cEjsi - 2c\4 s i ~ 17s 2 44 
+2s 2 44 + 47c?c|si + 14c?c|sis 2 + 14c? 4 + 36e? 4 

-Wc\4 s 2 + 10c?ctsi (38) 

As before, it can be seen that the denominator of the metric is given by the following trigono- 
metric polynomial: 

: = -21044 + 252c?c|sis 2 + 460e?cf - 1190c? 4 + 460c? cf 

-1190c? 4 ~ 45c i c 2 - 4 5c?C2 + 45c? 4° ~ 250c? 4° + 490c l c 2° 
-405c?c^° - 405c?°cf + 490c?°c 2 i + 45c?°c^ - 250c|°c 2 ! 



-584c?c£sis 2 + 120c?cisiS2 + 120c?c^sis 2 - 584c?c 2 , sis 2 
+ 10S1S2C?C2 + IOS1C2S2C1 - 344c?cfsiS2 - 120c?C2SlS2 
+332c?c|sis 2 - 344c?c^siS2 - 120sicf s 2 c? + 332sicls 2 c? 
+H80c?c| - 210c? 4 ~ c i° + 808c?cjsis 2 + 910c?c| 



-4° + 121c?°4° + 122c?c|sis 2 (39) 



We see that the determinant of the metric tensor remains non-zero in the space of power factors 

and thus defines a non-degenerate intrinsic geometry on the surface of fluctuation of phases. 

Finally, we may easily obtain the underlying scalar curvature which also has a rational 
polynomial form. It turns out that the scalar curvature can be reduced to the following form: 

„ (1)C . (2)C . (3)C . (4)C 

R= fiQ "V +"k +n R +n R U(]) 
n a c 1 r 2y2 (1)C . (2)C . (3)C V ' 

Our computation shows that the numerator of the scalar curvature takes the following trigono- 
metric polynomial expressions: 

n ( R )C := -495c? c| - 234c? c| -495c? 4 - 234c? 4 ~ 66c? 4° + 162c? 4° 
+13566c?4° - 91736c? 4° ~ 91736c?°cl + 13566c?°c^ - 66c?°c| 
+162c?°c 2 ' + 792c?c^sis 2 + 792c?c|sis 2 - 21064c? cfjsiS2 + 220c? c^s^ 



14 












V 



-8 -6 -4 -2 2 4 6 8 10 

a 

Legend 

Curve 1 



Figure 9: The determinant of the metric tensor plotted as the function of the equal power 
factor a :— a\ — a 2 , describing the complex power fluctuations in electrical networks. 
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+492cfc|sis 2 - 21064c? c£sis 2 + 220sie 2 's 2 c? + 492sic|s 2 cf 
+25506cf c| - c} 2 - c^ 2 + 2c^ 4 + 1152c 7 1 c 7 2 s 1 s 2 + 2c} 4 - 924c^c| 
+175360c} 1 c^ 1 s 1 s 2 + 233534c} 4° + 85296c? c%s lS2 + 186si c]; 3 c§ 
-1960sic} 3 c| + 5512sic} 3 c^ + 3214s!c} 3 cl - 34758sic} 3 4° 
+49084sic} 3 c^ 2 - 2054s 2 c?c 2 1 + 10386s 2 c?c^ - 9732s 2 c?4 1 
-33034s 2 c}°4 1 + 72538s 2 c} 2 c 2 1 - 38450s 2 c} 4 4 1 + 14436s 2 c? c\ 

n^ )C := 90s 2 c 2 c^ - 5796s 2 cfcJ + 594s 2 c?c| - 3878s 2 c}°c| 

+14424s 2 c}°c 2 ' - 6944s 2 c}°c 2 l + 6996s 2 c} 2 c| - 12292s 2 c} 2 c 2 ' 
-20348s 2 ci 2 cl - 3600s 2 c} 4 ci + 2428s 2 c} 4 c 2 ' + 18194s 2 c} 4 c 2 ) 
-14436sic?c^ + 6944sic?c 2 ° + 20348sic?c^ 2 - 18194sic?c^ 4 
-IOSSGsici 1 ^ + 9732sici 1 c| + 33034sic} 1 4° - 72538sic} 1 c 2 2 
+38450sicJ 1 c^ 4 - 5074s 2 c?c^ - 902s 2 ci 2 c| + 694s 2 c} 4 c 2 1 
+5074sic?c^ +440s 2 c 4 c| l - 30s 2 c} 4 c 2 + 2054sic} 1 c 2 
+5796sicjc| - 14424sic?4° - 2428sic?4 4 - 440sic?c 2 ' 
-660sic^c| + 12292sicjc^ 2 + 902sic?c^ 2 - 6996sic?4 2 

n^ )C := +3878sicfc 2 ° + 3600sic?c^ 4 - 694sic?c^ 4 + 242s 2 cJ°c| 
-594sic?c! + m0s 2 c\cl - 242sicf c\° - 20sicic^ 2 
+30sicic 2 4 + 20s 2 c} 2 c 2 - 90sic"cl - 49084s 2 c} 2 4 3 
+34758s 2 c}°4 3 + 1960s 2 c 4 c^ 3 - 3214s 2 c?4 3 - 5512s 2 c?4 3 
-186s 2 c 2 c 2 3 - 32s 1 s 2 c?c 2 1 - 288sis 2 c?4 3 - 5720sis 2 cf 
+5288sis 2 c?c 2 3 + 44672sis 2 c?c 2 :1 - 28816sis 2 c?c 2 3 
-132516sis 2 c?C2 1 + 72372si s 2 c\cf - 132516sis 2 c"c| 
-84592sis 2 c} 1 c 2 3 + 72372sis 2 c} 3 c 2 1 - 84592sis 2 cPc^ 
-288l6s 1 s 2 c\ 3 c 7 2 - 32s i s 2 c} 1 c 2 1 - 288sis 2 c} 3 c 3 ! 

n^ )C = -5720sis 2 cJ 1 c| + 44672sis 2 cJ 1 c5 + 5288sis 2 cfc| 

+12«is 2 c} 1 c 2 - 12sis 2 c} 3 c 2 - 12sis 2 cic 2 3 + 12s 1 s 2 cic 2 1 
+14504c?c^ 4 - 1552c? 4 4 - 26392c 2 2 c? + 122669c^ 2 cf 
-264682c^ 2 cJ° + 269177e 2 2 c} 2 - 102868c^ 2 cJ 4 - 26392c} 2 el 
+122669c} 2 c| ! - 264682c} 2 c\° - 1552c} 4 c\ + 14504c} 4 c| 1 
-56702c} 4 c| + 110054c} 4 c^° - 2sic} 3 + 6c} 4 ci + 6c?4 4 
+70c?c 2 2 + 1771c} 2 c 4 - + 70c} 2 c?. + 1771c 4 cf + 2s 2 cf 



-102868c} 2 c^ 4 - 56702c? cf + 110054c}°4 4 + 36560sis 2 c} 3 c 2 3 

-21276sic} 3 c^ 4 + 21276s 2 c} 4 c 2 3 + 36556c} 4 c 2 4 (41) 

We see further that the denominator of the scalar curvature has the following trigonometric 
polynomials: 

r^ )C := -210cfci + 252cf c|sis 2 - 68c 4 cf + 3660c? c§ 
-68c? c 2 + 3660c? c\ ~ 45c 2 c| - 45c? c\ 
-31c? cf + 991c? cf - 5074c? cf + 5257c? cf 
+5257c}°c| - 5074e[ c§ - 31c}°c| + 991c 10 c^ 
+172c?c^sis 2 + 120cjc|sis a + 120c? c^s^ + 172c?c|sis 2 
+10si« 2 c?c 2 + 10sic 2 s 2 ci + 5916c? c 2 sis 2 + 138c?c 2 sis 2 
-2132c? clsis 2 + 5916c? c^sis 2 + 138sic^s 2 c? - 2132sic;js 2 c? 
-10970c? c| - 210c? c\ - 672sic|c? + 672s 2 c|c? - c}° 
-4112c?c^is 2 + 70c? c\ - cf + 6836c} 1 c 2 1 sis 2 + 8265c}°c 2 ° 

r$° := -310c?c|sis 2 - 292s 2 c 4 c 2 1 + 2840s 2 c?c 2 1 - 7846s 2 c?c^ 
+8556s 2 c}°C2 1 - 3240s 2 c} 2 c 2 1 + 8876s 2 c?c| - 20s 2 c?c 2 1 
-220s 2 cfcl - 1596s 2 c?c| - 120s 2 c}°c| + 6560s 2 c}°c| 
-14604s 2 c}°c| + 1116s 2 c} 2 c| - 4784s 2 c} 2 c| + 6942s 2 c} 2 c 2 l 
-8876sic?cl + 14604sic?c^° - 6942sic?c 2 2 - 2840sic} 1 c| 
+7846sic} 1 cl - 8556810^02° + 3240sic} 1 c^ 2 - 688s 2 c?cl 
-24s 2 c} 2 cl + 688sic?c^ - 618s 2 c?c^ + 292sic} 1 c| 
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Figure 10: The curvature scalar plotted as a function of the power factors csi, 02, describing the 
complex power fluctuations in electrical networks. 



+220s lC l4 - 6560sic£4° + 618sicfc^ + 2080sic{c 2 i 
+4784sic5"4 2 + 2As 1 c{c\ 2 - 1116sicfc 2 2 + 120sicf4° 

:= -284s 2 cJ°c| + 1596sicfc| - 2080s 2 cfe 2 ' + 284sic?C2° 
+Ws 1 c 1 c\ 2 - 10s 2 ci 2 c 2 + 20s 1 c\ 1 c? 2 - 258cf cfsi 
-20ciC2°si + 258cfc 2 ! s 2 - 492c5'c 2 'si + 492cf cjs 2 
-196sis 2 ci 4 1 + 748sis 2 c5 i 4 1 + 768sis 2 4c 2 1 - 6886sis 2 c?C2 1 
-6886sis 2 ci 1 c| - ims 1 s 2 c\ 1 cl + 748sis 2 c} 1 c| + 768sis 2 eJ 1 c 2 ' 
+10sis 2 ci 1 c 2 + 10S1S2C14 1 + 7944 2 cf + 33664 2 e? 
-103034 2 cJ° + 6835c^ 2 c} 2 + 794c} 2 c 2 i + 3366c} 2 c| 
-10303c} 2 4° + 53c? e^ 2 - 489c{ 2 c% + 53c} 2 c| - 489cfc 2 2 

+2c^s 2 - 2c} 1 s 1 + 92c? c^ + 20cJ°c 2 s 2 - 92c? c|si (42) 



For the choices V = 1 and i?o = 1, the Fig. (JH]) shows that the determinant of the metric 
tensor has a commutative effect of the real and complex power flow fluctuations. This plot 
explicates the nature of the stability of a joint power flow in a realistic electrical network. 
The corresponding plot for the scalar curvature is depicted in the Fig. ljTO]) . This plot shows 
the global nature of a combined real and imaginary power flows in the electrical network. 
This analysis remains under the effect of Gaussian fluctuations of the realistic network of the 
resistance, reactance and impedance. 

For the equal network parameters a± = a and CL2 = a, the surface plots of the determinant 
of the metric tensor and scalar curvature are respectively shown in the Figs.flH]) and ffTTj) . We 
observe herewith that the stability of the electrical network power flow exists in certain distorted 
bands. These distortions are strong enough, so that they are capable to modulate the global 
properties of the fluctuation of the power flow. Specifically, we notice, for equal values of the 
network power factors, that the global instabilities exit for three specific values. 
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Figure 11: The curvature scalar plotted as a function of the equal power factor a 
describing the complex fluctuations in electrical networks. 



ai = a 2 , 



6 Conclusion and Remarks 

We have studied a power system planning, where the voltage stability is the main concern. In 
this setup, advancements in the electricity market, thus introduced, involve radical changes in 
the structure of the power systems. In the electrical environment of the power industry, we 
have considered an intrinsic geometric model to make a power system non-linearly efficient. 
Our model gives promising optimization criteria to select the optimal network parameters. 
The robustness of the model has been illustrated by variation(s) of the impedance angles, viz., 
phases. For a set of impedance, resistance and capacitance, our construction describes a definite 
stability character, with respect to the power fluctuations of the network(s). As a function of the 
trigonometric polynomials, we have demonstrated that the canonical fluctuations can precisely 
be depicted without any approximation. 

In the present paper, we have analyzed the statistical fluctuations in the real and imaginary 
power flows and thus characterized the network configurations. Furthermore, for the joint effect 
of the real and imaginary power equations, we have shown that the intrinsic geometric notion 
offers a clear picture of fluctuating network parameters. Such a configuration, as the limit of 
an ensemble of the power factor fluctuations, reduces to the specific electrical network. The 
present analysis does not stop here, but it invades the nature of the underlying imaginary power 
flow and explicates the global stability for the identical power factors. Our study thus offers an 
appropriate network design towards the stability of the existing linear optimization techniques. 

Finally, our proposition offers suitable tests towards the planning of the network parameters 
of finite component power systems. In our analysis, the intrinsic geometric model takes an 
account of the non-linear effects, arising in the stochastic power systems. A novel approach is 
herewith made possible in the history of network designs, and thus the present investigation 
provides appropriate network designs. The criteria deduced from our method can be used 
for determining optimization constraints both in the economic analysis of power systems, as 
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well as for setting up the operating point (s) of the power system, viz., network planning and 
compensation techniques. This task is left for the future investigation. 
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